function value = g(p)
x = p(:,1); y = p(:,2);
value = 1.0 ./ ((x-0.5).^2+(y-0.5).^2+0.01)...
        - 1.0 ./ ((x+0.5).^2+(y+0.5).^2+0.01);
end 